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Abstract 

We theoretically describe the quantum properties of a large Fresnel number degenerate optical parametric oscillator (DOPO) 
with spherical mirrors that is pumped by a Gaussian beam and is tuned at the subharmonic frequency to a given transverse 
mode family. We first analyze the classical problem and find that only the Laguerre-Gauss modes with lowest orbital angular 
momentum (QAM) are amplified above threshold. The transverse symmetry of the classically emitted signal field depends on 
the family index /: If / is even the lowest available QAM is zero and emission occurs in a radially symmetric mode; on the 
contrary, if / is odd the lowest available OAM is 1 and the emitted signal field has the shape of a hybrid Laguerre-Gauss mode 
(a linear combination of the two Laguerre-Gauss modes with OAM equal to ±1) that breaks the rotational invariance of the 
system. Next we focus on the squeezing properties of this DOPO model. As for the modes with lowest OAM we demonstrate 
that their quantum properties (in the linear approximation) are equal to the standard single mode DOPO (for even /) or to 
the recently analyzed DOPO tuned to its first transverse-mode family [Phys. Rev. Lett. 100, 203601 (2008)]. Concerning the 
rest of (classically empty) modes (having larger OAM) we find that combinations of Laguerre-Gauss modes with opposite OAM 
(hybrid Laguerre-Gauss modes) exhibit quadrature squeezing. This property is independent of the even or odd character of the 
family index and hence has nothing to do with the symmetry of the classically emitted signal field. Noticeably the amount of 
squeezing does not depend on the pump level (it is thus noncritical squeezing) and can be arbitrarily large for the lower OAM 
nonamplified modes. 

PACS numbers: 42.50.Dv, 42.50.Lc, 42.50.Tx, 42.65.Yj 



I. INTRODUCTION 

Degenerate optical parametric oscillators (DOPOs) are 
nowadays the standard squeezed light source. Let us re- 
mind that a light mode is said to be squeezed if the fluc- 
tuations in one of its quadratures are below the stan- 
dard quantum limit, which is defined as the vacuum 
fluctuations level of that quadrature [J, [2) DO- 
POs, squeezing is accomplished thanks to the parametric 
down-conversion process occurring in the nonlinear crys- 
tal together with the interference between the intracavity 
fleld and the external vacuum fluctuations that enter into 
the cavity through the output mirror . Quantum noise 
reductions as large as lOdB (90%) have been experimen- 
tally demonstrated with DOPOs 

High precision measurements areperhaps the best 
known applications of squeezed light '311 , but applications 
to quantum information with continuous variables are 
becoming increasingly important [61 , as squeezed light is 
the essential ingredient in generating continuous-variable 
entanglement. Improving the quality and reliability of 
squeezing is thus an important goal, but the genera- 
tion of squeezed light with particular spatial distributions 
could also be important. This has been shown to be of 
utility in, e.g., high precision positioning [7|. Recently 
Laguerre-Gauss beams are attracting much attention be- 
cause of the many potential applications of the orbital 
angular momentum (OAM) carried by these light beams 
Q , and the generation of non classical Laguerre-Gauss 
modes could thus lead to new phenomena in the interac- 
tion between these light fields and matter. 

Here we consider squeezing generation by means of 
DOPOs with large Fresnel number cavities composed of 
spherical mirrors. In such systems, the cavity can sustain 



the nonlinear interaction for several transverse modes 
(diffraction losses are ideally suppressed) giving rise to 
new results concerning squeezing. Quantum fluctuations 
in large Fresnel number cavities have been studied in 
the past, and new phenomena resulting from the inter- 
play between quantum fluctuations and transverse pat- 
tern formation have been predicted. Concerning cavities 
composed of planar mirrors, the phenomena of quantum 
images ^] (below threshold) and the perfect noncritical 
squeezin g of the emerging transverse pattern linear mo- 
mentum [l^, [ll| (above threshold) have been predicted. 
Regarding spherical mirrors, only below threshold op- 
eration has been studied [l^, showing that also in this 
case there appear quantum images anticipating the above 
threshold pattern. 

Here we consider a DOPO cavity with spherical mir- 
rors that is pumped by a Gaussian mode. The novelty of 
our work with respect to previous analyses of large Fres- 
nel number DOPOs consists in that we study the sys- 
tem above its oscillation threshold, while previous works 
have treated the below threshold case j^, ■ This has 
lead us to quantum fluctuations phenomena that have 
not been previously described. Specifically, we shall as- 
sume that the DOPO cavity is exactly tuned to a partic- 
ular family of Laguerre-Gauss transverse modes labeled 
by an integer / (consisting of / -I- 1 frequency-degenerate 
Laguerre-Gauss modes with OAMs ±/, ± (/ — 2),...,±Zo 
with = or 1 for even or odd / respectively; this is 
reviewed in Appendix A). After deriving in Section II 
the system's model, we first demonstrate (Section HI) 
that the DOPO above threshold emits a signal field con- 
sisting of pairs of photons with opposite OAM, -f/g and 
— ^o; i.e., the DOPO emits in the transverse mode with 
the lower possible OAM, a result not derived previously. 



The squeezing properties of this mode turn out to coin- 
cide with those derived in the past, as the case Iq = 
is equivalent to the standard DOPO model 0, and the 
casejn = 1 generahzes what has recently been described 
m [ilQ. After reviewing the squeezing properties of 
cases /o = and 1 in Subsection IV. A, we concentrate in 
the squeezing properties of the rest of modes (those with 
I > Iq that remain off when the DOPO is above thresh- 
old) in Subsection IV. B. In this last subsection we will 
demonstrate the main result of the present work, namely 
that certain combinations of opposite 0AM modes ex- 
hibit large squeezing (which is larger for the larger val- 
ues of / and the lower values of /). Remarkably, this 
squeezing is noncritical, i.e., is independent of the pump- 
ing value and admits a simple explanation (Subsection 
IV. C). We find that this new result is relevant because 
it establishes a simple way for generating squeezed vacua 
with shapes different from the Gaussian or the TEMio 
modes. In Section V we review the main conclusions of 
our work. 



II. DOPO'S QUANTUM MODEL 

We consider a type I DOPO with spherical mirrors 
pumped by a coherent Gaussian beam of frequency 2ujq, 
matched to the fundamental transverse mode at that 
frequency. Within the cavity, pump photons are down- 
converted into signal photons of frequency coo at a x^^"* 
crystal that is placed at the cavity waist. The impor- 
tant assumption of this work is that the resonator is 
tuned so that luq coincides with the resonance frequency 
of a unique transverse mode family /, and is detuned far 
enough from any other family. Within these conditions, 
the total electric field at the resonator waist plane can be 
written as 

E{r,t)=Ep{r,t)+Es{r,t), (la) 
Ep (r, t) = iTpAp (r, t) e-^'-^"* -I- H.c, (lb) 
4 (r, t) = iJ^sAs (r, t) e-'^"* + H.c, (Ic) 

where we have introduced the single-photon field voltages 



eoHcLcs ' 



(2) 



being is the crystal refractive index and Lcs the effec- 
tive resonator length (see Appendix A), and the slowly 
varying envelopes 



Ap (r, t) = aoo (t) (r) , 



A. 



1 



(l + <5o/) 



a+i 



(3a) 
(3b) 



where the interaction picture boson operators obey the 
canonical commutation relations [a™ (<) , a}^ (<)] = Smm 



the Kronecker symbol Sqi is introduced because a_o {t) = 
a+o (t), and ^'^'_;^y2 resonator Laguerre- 

Gauss modes at the resonator waist plane. The expres- 
sions for these modes as well as some of their properties 
are given in Appendix A. Here we just note that 



(4) 



where r = (r cos 0, r sin 0) are the resonator waist plane 
cartesian coordinates (r and are the corresponding po- 
lar coordinates), A/"^ is a normalization factor. 
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(5) 



LJj is the modified Laguerre polynomial, and w is the 



beam waist radius (that depends on frequency). In (j3bp 
the sum in I runs along all the members of the transverse 
mode family f [I G {/, / — 2, Zq}] and this will hold 
for all sums in I along the rest of the article. Finally, 
note that the Laguerre-Gauss modes appearing in (j3ap 
and (j3bp are evaluated at pump and signal frequencies, 
respectively, which determine the beam waist radius w 
(see Appendix A). 

For the sake of later use we mention here that instead of 
the Laguerre-Gauss modes one can use the hybrid modes, 
defined (for I ^ 0) a,s 

(r) = -i= (r) + (r)] = V2Mlu\ (r) cos {1^) , 

(6) 

(r) = -i, (r) - vj/;' (r)] = V2AA>i (r) sin (/0) , 

(7) 

(see Appendix A) in order to expand the slowly varying 
envelope of the signal field, whose corresponding boson 
operators relate to the Laguerre-Gauss mode ones as 



acd = [a+i (<) + a-i {t)] /v2, 
fls,; = i [a+i (t) - a^i {t)]/V2. 



(8b) 



For / = we still use the Laguerre-Gauss mode 5'° (r). 

In the interaction picture, assuming perfect phase- 
matching as well as exact resonance between the fields' 
frequencies and the cavity resonances, the Hamiltonian 
of the system is _ff = iJoxt + -ffint , with 



p {^\) 



00 - ooo ) , 



(9a) 



sr^ XI 



^1 



30/ 



(^aooa^jfll; - aJoO+ia.;^ , (9b) 



where J¥ext describes the external pumping process and 
ifint describes the down-conversion process occurring in 
the x*^^-' crystal. Notice that 0AM conservation imposes 
the creation/annihilation of a pair of signal photons each 
with an opposite value of /. 
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The nonlinear coupling constants xi in ^int read 



XI = 12 



(2)7 

It. 



UJQ 



JlrL. 



cff 



3/2 



Treo 



(10) 



with x*-^' the second order susceptibility of the nonlinear 
crystal whose thickness is Ic, Wp the beam spot size at 
the pump frequency, and 



j as Xj = (a^ + a^ ) and Yj 



{at 



aj) . Outside the 



cavity, it is measured the variance spectrum of quadra- 
ture Xj (analogously for Yj), which can be calculated 
as V°j^ (w) = 1 + 5°"* (w), with the squeezing spectrum 
given, as a function of the intracavity amplitudes, by Q 



27, 



+ 00 



dre 



{X,{t),X,{t + T)), (14) 



[(/-0/2]! 
[(/ + 0/2]! 



+ 00 



due 



L 



(/-0/2 



(U) , (11) 



which are proportional to the overlapping integrals be- 
tween the three modes involved in the particular para- 
metric process. This means that the nonlinear coupling 
between pump and signal photons is larger the lower is 
the 0AM of the latter ones, i.e., xj < Xj-2 < ■■■ < Xio- 
This property will play an important role as we show be- 
low. Let us finally remark the the pump parameter £p 
is proportional to the external pump amplitude and is 
taken as real without loss of generality. 

We will be interested in calculating normally ordered 
correlations of different mode operators, and to do so we 
use the generalized P representation and its equivalent 
set of Langevin equations [lB| . In this representation to 

every pair of boson operators ^a, , a^^ it corresponds a 
pair of independent stochastic amplitudes (aj,Q;^) that 
are complex-conjugated in average, i.e., (ctj^) = {o-j)* ■ 
The stochastic average of any function of (a,, a^) equals 
the corresponding normally ordered, quantum mechani- 
cal expected value. The equations of evolution of these 
amplitudes are derived in Appendix B by following the 
standard procedure [l^, [l3|. Assuming that losses oc- 
cur just at one of the cavity mirrors at rates 7p for 
the pump mode and 7s for all possible transverse signal 
modes (hence the assumption of a large Fresnel number 
resonator) the Langevin equations of the system read 



c sr^ Xi 

aoo = Op- jpUoQ - > — -^a+ia_i. 



an 



00 ^ Ip^oo' 



XI + + 



Ooi 



(12a) 



(12b) 



+ 



a 



I ^ ~'^sa±i + Xia'^i^oo + ^/xiaoo£.±i {t) , (12c) 
I = -Isa^i + Xia^i^oo + \/xiaoo^±i (*) , (12d) 



where Z e {/, / — 2, ^o}- The noises are inde- 

pendent complex Gaussian white noise sources verifying 

(6 (t)) - (t)) = and 

(6 (t) (^')> = (t) fct it')]*) = 5u'6{t t') , (13) 

being null the rest of correlations. The rest of noises 
verify U (t) = ^ C*) and e+ (t) = [C+ (t)] * 

Within the generalized P-representation, we define the 
amplitude and phase quadratures, respectively, of a mode 



where the factor 27, comes from the input-output rela- 
tions [ll|. We use the notation (a, b) = {ah) — (a) (6). No- 
tice that 5*°"* (ws) = — 1 signals perfect squeezing outside 
the cavity for quadrature Xj at detection frequency ujg 
(which must not be confused with the optical frequency, 
as it has contributions of every pair of modes lying in 
opposite sidebands around the optical frequency ujj +0Js, 
where uoj is the carrier frequency of the detected mode 
[4j). Hence, by solving the Langevin equations for the 
amplitudes aj in terms of the noises ^fe, one is able to 
obtain the variances of the quadratures involved in the 
problem, and thus the squeezing properties of the field. 



III. CLASSICAL EMISSION 

The classical equations describing the field inside the 
DOPO are obtained from the quantum Langevin ones by 
identifying the stochastic amplitudes aj with the classical 
normal variables of each mode, by making a'j — > a*, and 
by neglecting noise terms. They read 

Qoo = - 7pCtoo - -, ^\ a+iot-u (15a) 



a±i = -lsOi±i + xioh^i^oo- 



(15b) 



Equations (jl5[) have two types of stationary solutions. 
If we define the normalized pump parameter as 



Xlo 

Ipls 



pi 



(16) 



it is easy to prove that the below-threshold solution 

Oioo = £php, a; = V;, (17) 

is stable for cr < 1 and unstable for ct > 1 (a = 1 thus 
defines the classical threshold for emission). Apart from 
this trivial solution, there are (/ — Iq) /2 + 1 possible sta- 
tionary solutions in which the signal field is nonzero. The 
form of these solutions is 



aoo ="/s/Xk, a±i=0 y I ^ k, 



a±k = Pke 



with 



Xlo 



2 
Pk 



K-k 



l±So£ 

g^Kk 



Xictoo 

Is 



Xk _ 

Xlo Ilo ' 



(18a) 
(18b) 

(19) 
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Figure 1.- Transverse profile of tlie signal field modulus 
above threshold. 



see Eqs. (jlOp and (fTT|) . Note that in any of these so- 
lutions only two opposite 0AM values are excited (±fc), 
the rest remaining below threshold. As for the phases 
Ok, Oq — and 0fe^o = ^ is arbitrary. This arbi- 
trariness appears because Eqs. p5|) have the symmetry 
a±i —^ a±i exp (±i/3) , reflecting the rotational invariance 
of the system, which leaves undefined the phase difference 
between opposite 0AM modes. 

As for the stability, it is easy to show that the only 
solution that is linearly stable above threshold is precisely 
that in which the lower 0AM modes are switched on, i.e., 
Eq. (fT8|) for k — Iq. These arc, in fact, the solutions with 
the lowest threshold a = 1 [see Eq. (US])], as xia > Xi 
already stressed [see Eq. ([T(I| ]. Plugging Eq. PS)1 into 
Eq. (j3bp one finds the slowly varying amplitudes of these 
stable classical solutions as 

As (r) = po*^/2 (r) for even / (20a) 
As (r) = N/2pii/i (r, q}-9) for odd / (20b) 

with pj given by Eq. ^5)1 . In Fig. 1 the square modulus 
of these slowly varying envelopes are shown for the first 
four families (note that we have chosen 9 = tt/A for these 
figures). 

A most important difference occurs between the two 
cases corresponding to / odd or even: For even / the clas- 
sical emission As (r) is rotationally symmetric as {^) 
is [see Eq. Q], while for odd / the emitted pattern 
breaks that symmetry. 



IV. SQUEEZING 

Now we perform the analysis of quantum fiuctuations 
affecting the classical emission above threshold. First, 
in order to make the calculations as simple as possible, 
we consider the limit jp ^ 7^ in which the pump mode 
can be adiabatically eliminated in Eqs. as aoo = 

^p/lp ~ 'i+Sn: ■ The resulting equations read 

a±i - Is {-a±i + mNa+i) + ^/^y^£.±i (t) , (21a) 

"±i = 7. i~a±i + i^iN+a+i) + /WV+C^ (t) , 

(21b) 

where 

N^a- g^Y^-^l—a+i.a.v , (22) 

^l + (30,i' 

and iV+ is as TV but replacing av by af,. Moreover, as 
we will study fluctuations in the linear approximation, 
the only nonlinearities contained in TV that will give a 
contribution correspond to the classically excited modes 
il' = Iq) and hence 

^---9'^^;^, (23) 

which is the expression we will use in the following. 

We are going to distinguish between the special cases 
I = Iq on one hand and cases with I > Iq, as their prop- 
erties turn out to be quite different. 

A. Special cases / = lo 

Here we analyze the squeezing properties of the modes 
that are classically excited. In this case ki = k/j, = 1 
[see Eq. P^ ]. and the Langevin equations (PTjl and (l^^)) 
become 

• Case lo — 

do = 7s {-ao + Na+) + y/^^o (t) , (24a) 
N = a- y^aoao, (24b) 

• Case Zo — 1 

a±i=ls{-a±i+Nal^) + y/^^±i{t), (25a) 
N = a- g^a+ia^i (25b) 

We note that in neither case the information about the 
family order / appears (up to here, that information ap- 
peared through the reduced coupling parameters ki [see 
Eqs. (HH), (Uni) and (HI])]). Moreover, Eqs. ^ coincide 
with those for the paradigmatic single-mode DOPO (see, 
e.g., 0), while Eqs. (|^5)) do with those for the recently 
studied DOPO tuned to its first transverse mode family 
[isl. [igj. The conclusions are straightforward: 
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(1) For any even / (Iq = 0) the quantum properties of 
the classically excited mode coincide with those of 
the usual single-mode DOPO in the linear approx- 
imation. At threshold (cr = 1) ideal squeezing is 
predicted in this approximation and this squeezing 
degrades as we move apart from threshold. That is 
why we speak about critical squeezing in this case, 
as the system parameters need to be tuned to a 
particular value for obtaining the optimum squeez- 
ing level. In mathematical terms, this is a conse- 
quence of the existence, at the bifurcation point, 
of a null eigenvalue and of a companion negative 
eigenvalue that reaches its minimum possible value. 
The former is the responsible for the " infinite" fluc- 
tuations of the amplitude quadrature of the signal 
mode (which is proportional to the eigenvector as- 
sociated to the null eigenvalue), while the latter is 
the responsible of the complete suppression (in the 
linear approximation) of the fluctuations in the sig- 
nal mode phase quadrature. (See below for a more 
comprehensive explanation.) 

(2) For any odd / {Iq ~ 1) the quantum properties 
of the classically excited mode coincide again with 
those of the single-mode DOPO, i.e., it shows per- 
fect squeezing only when working exactly at the 
bifurcation. Unlike the previous case the classi- 
cal emission breaks the rotational invariance of the 
system and then, apart from the previous criti- 
cal squeezing property, the system exhibits perfect 
squeezing in the remaining mode that is spatially 
crossed with respect to the classically emitted one 
[T3 |. As the reason for this lays on the rotational 
symmetry breaking (and then on 0AM conserva- 
tion (isjl. it has nothing to do with bifurcations and 
the predicted squeezing is noncritical: The same 
squeezing level (perfect, in this case) is obtained 
at any pumping level. From a mathematical view- 
point one of the eigenvalues governing the stability 
of the classical solution above threshold is always 
zero (the eigenvector associated to this eigenvalue 
is said to be a Goldstone mode), reflecting the inde- 
terminacy of the phase difference between the two 
signal modes [see Eq. (ITSl) ]. This phase difference 
is nothing but the orientation of the emitted hy- 
brid mode (the TEMio mode in [l^) that results 
from the coherent superposition of the ^'^"'_i)/2 ^^'^ 
^(/'-i)/2 Thus the physical meaning of the 
Goldstone mode is that the orientation of the signal 
mode emitted above threshold by this DOPO dif- 
fuses with time, and is thus undetermined. This 
eigenvalue is accompanied by another eigenvalue 
which always takes its minimum possible value, 
with the consequence that the amplitude of its asso- 
ciated eigenvector has no fluctuations at all. This 
last eigenvector can be easily identified with the 
OAM of the emitted pattern, which is nothing but 
another hybrid mode spatially crossed (orthogonal) 



with respect to the bright one. This shows that 
the generation of noncritical squeezing through the 
spontaneous rotational symmetry breaking mecha- 
nism is not particular of the model considered in 
[isj , but should be a ^general phenomenon. We ad- 
dress the reader to [l3] for full details of this phe- 
nomenon. 



B. Cases I > lo: Squeezing of higher OAM modes 

In this case the amplitudes ai and a'^ are very small 
(they are null in the classical limit even if the system is 
above threshold) and then, in the linear approximation 
we are considering, the nonlinear function N (|23p appear- 
ing in (PT|) must be evaluated at its classical value as 
not doing would introduce higher order corrections that 
must be neglected in the used linear approximation. One 
then obtains A^ = 1 and the Langevin equations become 

a±i = Is {-a±i + Kia'^i) + Vtv«7^±/ (<) , (26a) 

= Is (-a±i + + Vl^^ti (*) ■ (26b) 

As is evident, the evolution of the different pairs of OAM 
±/ is decoupled, what allows us to analyze separately the 
quantum properties of each I couple. Moreover, the evo- 
lution of these fluctuations does not depend on whether 
the family is odd or even. 

The study of Eqs. (|26|) is facilitated by expressing 
them in vector form as 



a; = ClSii + VTsKzI; (t) , 



with 



a; col a-i, , 



col 



(6 (o,e,+ (o,erw,fc+w]' 

/-I 

-1 Ki 
Ki -1 

\ki -h 



(27) 

(28a) 
(28b) 

(28c) 



where we remind that < < 1 [see Eqs. (|19p ]. 
Linear matrix Ci has eigenvalues X\ 



A 



-7s (1 - Ki) and = A^i 
sponding eigenvectors [IT 



-7s {1 + Ki) with corre- 



4,2 = ^C0l(l,±l,±l,l), 



[5,4 = -col(l,±l,Tl,-l). 



(29a) 
(29b) 



By projecting the linear Langevin system (|27p onto 
we find 



^2,3 



~7s (1 T Ci + \/7lKi^i,4 W ' (30a) 
-7s(lTKi)c3 + «V7^'?2,3(0: (30b) 



5 



where we defined the projections 

4 {t) = w^. • ai (t) , 
and the four new real noises 



(31) 



4a (t) = • a (t) - Re [^i (t) ± e+ it)] , (32a) 
V2,3 it) = -«w^,3 • (t) = Im [C, (t) T C (*)] - (32b) 
which satisfy the statistical properties 

{V',{t))=0, {vLitl)vUt2))^SrnnSih-h). (33) 

Equations ([50)1 are solved with standard methods. The 
variance spectra for projections Cj, defined as 



+00 



dre 



{c,it),C,it + T)), 



read, in the stationary limit t 3> A ■ 



(34) 



(35) 



The advantage of the eigensystem method we are us- 
ing [l^, [ll|, [1^ is that the relevant quadratures of the 
problem appear in a natural way, as the projections Cj 
are related to the problem quadratures through 



•2 7 



where 



V2c 



3' 



(36a) 
(36b) 

(37) 
(38) 



refer to the hybrid mode ^/j (f-i)/2 ('") '"^^^ easily 
seen from Eqs. ([S]). These simple relations are the con- 
sequence of the appropriate choice of the eigenvectors of 

We are now in conditions to find out the squeezing 
spectra of the quadratures. By using Eq. to- 
gether with the relations above, we get that 5°"*; [uo) — 

SZM = ^lsC\{^) and S--1,{lo) = 5-^(0;) = 
4756*3 (w), i.e., the squeezing of the phase quadratures 
are 



5°"';(c^) = 2, j = c,s. 



(39) 



This expression clearly shows that the phase quadratures 
y^.i U = c,s) of both i?^^(^_,)/2(r) and i?^_(^_,)/2 (r) 
hybrid modes (i) have the same fluctuations properties, 
(ii) are noncritically squeezed as S°^ji (uj) < at any 
noise frequency uj and this value is independent of the 
distance from threshold [it is however dependent on the 
0AM value through the value of the ratio ki (fT5|) ]. (iii) 
and exhibit maximum squeezing at w = 0. This is a main 
results of our article and we pass to discuss it. 




Figure 2.- Variance spectra of the phase quadratures of some 
Hybrid modes with I ^ lo- The inset shows the same 
quantity at a; = as a function of kj. 



C. Discussion 

We have shown that the light emitted by a DOPO with 
a large-Fresnel-number cavity with spherical mirrors ex- 
hibits a variety of squeezing properties. On the one hand 
there is the usual critical squeezing appearing at the os- 
cillation threshold associated with the classically emit- 
ted solution. There is also the noncritical squeezing due 
to the spontaneous breaking of the rotational symmetry 
that appears above threshold for odd / (^o = 1) [13j . 
Also a new type of squeezing behavior has been derived 
for all hybrid modes with I > Iq, which remain off at the 
classical level once the oscillation threshold is crossed, 
consisting in that their phase quadratures are noncriti- 
cally squeezed for any I, as shown by Eq. ([55)1 . 

The variance spectrum for either of the phase quadra- 
tures V°f (u) = 1 + S°'^li (uj) for j = c or j = s, is 
represented in Fig. 2 for different values of ki [see Eq. 
pt| ]. In the figure, V°f (uj = 0) is also shown, at the 
inset, as a function of k;. Clearly, as k/ = Xi/Xh ~^ 
what occurs as ^ ^ Zq, nearly perfect squeezing is found 
at u; = 0. Notice that more than 90% of squeezing can 
be achieved if k/ > 0.5. But the ratio ki = Ii/Iig, which 
is obtained from Eq. PH)) . is solely determined by ge- 
ometrical reasons, and thus the squeezing properties of 
these higher 0AM modes do not depend on the system 
parameters (of course we must remind that we assumed 
equal cavity losses for all modes). 

In order to better appreciate the large amounts of 
squeezing exhibited by these modes, we give their as- 
sociated noise reduction in percentage for different even 
(odd) families in Table 1 (2). Notice that the largest 
squeezing occurs for large values of / and small values of 
Z, and works better for even families. 

It is to be remarked that both the hybrid mode 7?' „ (r) 
and its orthogonal „ (r) have the same squeezing prop- 
erties. This means that the orientation of the mode is 
irrelevant, as an hybrid mode rotated an arbitrary angle 



6 



Table 1 




f\l 


2 


4 


6 


8 


10 


12 


1 


88.9 












4 


96.0 


49.0 










6 


98.0 


71.0 


18.1 








8 


98.8 


81.6 


36.8 


5.6 






10 


99.2 


87.4 


51.4 


14.7 


1.6 




12 


99.4 


90.9 


62.1 


24.9 


5.1 


0.4 



Table 1.- Percentage of noise reduction for the non amplified 
Hybrid modes I 7^ laying in even families. Note that large 
levels of squeezing are obtained for the lower angular 
momentum modes. 



Table 2 




f\l 


3 


5 


7 


9 


11 


13 


3 


75.0 












5 


88.9 


33.1 










7 


93.8 


55.6 


10.8 








9 


96.0 


69.1 


24.9 


3.1 






11 


91.2 


77.6 


38.0 


9.1 


0.9 




13 


98.0 


83.0 


49.0 


16.6 


3.0 


0.2 



Table 2.- Same as Table 1, but for odd families. Large levels 
of squeezing are also obtained in this case for the lower 
angular momentum modes. 



P respect to the x axis, which is given by 

HlA^)^Hl^{r,<l,~p) (40) 
= Hl^ (r) cos (W) + Hl„ (r) sin (1(3) , (41) 

also has the same squeezing properties. This is in clear 
contrast to the perfectly squeezed mode in the I = Iq = I 
case [l3| , which has not an arbitrary orientation but is or- 
thogonal to the classically excited mode at every instant. 
Of course, the reason is that the modes here analyzed are 
below their threshold, at difference with the case in p^ . 

So far we have proven that the higher 0AM modes 
show noise reduction in their phase quadrature. Let us 
see why is this occurring. The coefficient that rules the 
squeezing properties of the I ^ Iq modes is k;, see Eq. 
((39)) . Both above and below threshold, k/ can be written 
in terms of the stationary pump field as k; = Xic^oo/ls, 
see Below threshold aoo increases linearly with the 
external pump amplitude [see Eq. (jl7p ] and hence so does 
Ki. But once the bifurcation is reached and the DOPO 
starts emitting the mode(s), the value of aoo saturates, 
and remains constant irrespective of the pump value [see 
Eq. HH])] leading to k; — Ii/Iiq above threshold, which 
is independent of the pump value. This implies that the 
value of K/ will remain fixed to that at the bifurcation 
once the DOPO is above threshold and, consequently, 
the squeezing level of these modes becomes noncritical. 



V. CONCLUSIONS 

We have shown that tuning large Fresnel number DO- 
POs with spherical mirrors to transverse families is a sim- 
ple way for generating squeezed light with the shape of 
hybrid Laguerre-Gauss modes with 0AM I > Iq, with 
Zo = or 1 depending on the even or odd character of 
the selected transverse mode family. The behavior of the 
system can be resumed as follows: Above threshold, only 
the hybrid mode with I = Iq is amplified, and the rest 
of modes (those with I > Iq) remain off. The amplified 
modes exhibit squeezing properties that^can be found in 
the literature (see for = and 13] for Iq = 1) as 
shown in Section IV A. 

The surprising result is that hybrid modes with I > Iq 
exhibit a large degree of squeezing together with the fact 
that this squeezing is noncritical, i.e., it remains fixed 
irrespective of the pump level when the system is above 
threshold. Thus, squeezed vacua with the shape of higher 
order hybrid modes can be generated by normal DOPOs 
if these have a large Fresnel number. 

In our model we assumed equal cavity losses for all 
signal modes (in fact, equal cavity losses for all signal 
modes belonging to the same family; this is the precise 
meaning we give to the expression large Fresnel num- 
ber cavity). As we have explained, the squeezing level 
of modes with 0AM I > Iq depends on the distance be- 
tween their threshold and that for the mode(s) with low- 
est 0AM Iq. Hence, differences in the values of their 
decay rates could modify quantitatively the results we 
have presented, as this would change the threshold for 
the different modes. If fact, if tailoring the cavity losses 
for the different modes would be possible, one could ob- 
tain even smaller quantum fluctuations for these empty 
modes if their corresponding thresholds were made closer 
to that of the mode. We also assumed perfect cavity 
resonance for a particular family of modes. We do not 
expect that detuning changes our main conclusions, as 
detuning affects equally all modes within the same fam- 
ily (apart, of course, from the fact that if the detuning 
is half the value the transverse free spectral range, com- 
petition phenomena between the involved families will 
manifest, something that is not considered in our model). 
Adding detuning would have the same consequences as 
in the single-mode DOPO: It will increase the threshold 
and change the phase of the amplified mode (and hence 
of the squeezed quadratures). 

We finally note that Laguerre-Gauss modes are be- 
coming important for several purposes, e.g. in quantum 
information processing and in manipulation of atoms ^ 
just to mention a couple. We hope that our results can be 
of relevance as squeezed hybrid modes or, equivalently, 
entangled Laguerre-Gauss modes, can be of utility for 
these purposes. 

This work has been supported by the Spanish Minis- 
terio de Eduacion y Ciencia and the European Union 
FEDER through Projects FIS2005-07931-C03-01 and 
FIS2008-06024-C03-01. C.N.-B. is a grant holder of the 
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(Spain) . 



APPENDIX A: CAVITY MODES 

For the sake of clarity, we find it convenient to review 
the main properties of the cavity modes in a Fabry-Perot 
resonator with spherical mirrors (see, e.g., [2l| for more 
details). Within the paraxial approximation, it is well 
known that the Laguerre-Gauss modes form a complete 
set of spatial modes describing the light inside the res- 
onator. Let Ri and R2 denote the curvature radius of the 
cavity mirrors, and L^s = L — (l — 1/rtc) Ic the effective 
cavity length, being L the geometrical length of the res- 
onator, and Ic and Uc the length and refractive index, re- 
spectively, of the x*'^^ crystal. Then the Laguerre-Gauss 
modes at the resonator waist plane can be written as 



ir)=K< (r)exp(±z;0)^ 



(Al) 



with r — (x, y) the transverse coordinates being r — 
r (cos (j), sin 0) its polar decomposition, A/'/j a normaliza- 
tion factor and 




g 

Figure 3.- Resonance frequencies of the first 4 families 
corresponding to three consecutive longitudinal modes, as a 
function of the g parameter of a symmetrical resonator. 
Different families corresponding to the same longitudinal 
mode have the same color and dashing. The family order / 
increases with the frequency. 



and the transverse part of the resonance frequency is 
given by 



(r) = - 

w 



1 / V2r 



w 



2r2 



(A2) 



Awn; = - — (1 + 2n -I- /) arccos (1/5152) . (A7) 

^cff 



being the modified Laguerre polynomial with radial 
and polar indices n, Z e N, which are given by Rodrigues 
formula 



n\ v' av 

By choosing the normalization factor as 



TT {n + /)! 



(A3) 



(A4) 



the following orthogonality relation holds 

/ d0 / rdr [vI/±' (r)] * (r) = 6w6r^n' ■ (A5) 

Jo JQ 

The beam spot size at the cavity waist, w, is given by 



2cLeff Vffig2 (1 -gig2) 
gi+92- 25152 



(A6) 



with 5j = 1 — Lcs/Rj and lu the beam frequency. 

The Laguerre-Gauss basis is recommended in order to 
visualize the QAM of the field, as these modes are eigen- 
states of the 0AM operator —id^ with eigenvalues ±1. 
Concerning the resonance frequency of the different 
modes, they are different in general for each mode. Con- 
cretely ujqni = <l''^c/ + Aujni, with q an integer (differ- 
ent q's correspond to different longitudinal cavity modes) 



Hence, cavity modes having the same family order 
f = (2n + 1) have the same frequency and are said to 
be members of the same family /. ft is clear that fam- 
ily / consists of the set of / -I- 1 Laguerre-Gauss modes 

{*f/-0/2}' "^^^^^ I = fJ- 2,-,^o, having 0AM ±/, 
± (/ — 2),...,±^o, respectively. The lower 0AM modes 
have Zo = or 1 for even or odd / respectively. In 
Fig. 3 we represent Wqni (in units of the free spectral 
range, fsr = ttc/Lcs) for several modes as a function of 
5 for symmetric resonators (those for which 51 = 52 = 5) 
which are stable (i.e., — 1 < 5 < 1; be reminded that 
5 = 1 corresponds to a planar mirrors resonator, 5 = 
to a confocal resonator, and 5 = —1 to a concentric 
resonator). It can be clearly appreciated that different 
transverse families corresponding to different longitudi- 
nal modes can have the same frequency for rational values 
of (arccos5) /tt. Below we assume that different families 
have different frequencies, i.e., that this quantity has an 
irrational value. 



APPENDIX B: LANGEVIN EQUATIONS 



Next we outline the derivation of the Langevin Eqs. 
(fT^ . Neglecting the effect of thermal photons, the mas- 
ter equation ruling the evolution of the system's density 
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operator [l6| reads 



+ ^^Is (^'^aigaj - ajaig - gajai 
±1 



(Bl) 



where losses are assumed to occur in only one of the cav- 
ity mirrors and H is given by Eg. ([9]). Now we use the 
positive P representation [1^, |T3| of the density opera- 
tor, which allows the evaluation of expected values of nor- 
mally ordered operators as averages of functions in phase 
space. By using standard methods [isl. Iv^. the master 
equation (jBip is then transformed into a Fokker- Planck 
type equation for P (a, q:+), where all the coherent am- 
plitudes are collected into the vector 



a = col (aoo, "oo^ "1;) 
The equation for P reads 



(B2) 



d_ 

dt 



P = 



2J 



P, (B3) 



where ai means the i-th. element of ct, with 

XI 



I 



V 



35 



A3 = -7sa+; + X/aooa|/, 



(B4) 

(B5) 
(B6) 
(B7) 



and A2^Afi are obtained from ^1,3 5 by changing aj < — > 
. Any other element of the diffusion matrix T) is null. 



Next, as stated by the Ito theorem [24I, this Fokker- 
Planck equation is mapped onto the set of coupled 
stochastic (Langevin) equations 



a = A{a.)+B{a.) ■ f] (t) , 



(B8) 



where 77 is a vector with real white noises as components, 
each of them satisfying the statistical properties defined 
in Eq. (|33p. A is a vector with components Ai defined 
above; and the noise matrix B is defined hy V = B ■ B^. 
The equivalence between the Fokker-Planck equation and 
the Langevin system has to be understood as (/ (a))p — 
(/(«)) 



stochastic ' 



i.e., phase space averages are equal to 
averages made by using the statistical properties of the 
noise vector t] {t) . 

As for the diffusion matrix 2?, it is composed of un- 
coupled minors for each one of the ±1 and 00 subspaces, 
and then so does matrix B. In addition, the noise matrix 
minor associated to the pump mode is the zero matrix as 
this is its associated minor in V. One possible choice for 
the minors associated to ±1 modes is 



(B9) 




for / 7^ and 



So = Vxa 




(BIO) 



for / = 0. Finally, inserting this noise matrix into the 
equivalent Langevin system (|B8|) . we get Eqs. p2|l . 
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